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ABSSUCT 


jpmMfBt investigations oosprises of an iBplementati<»i 
of a oorapnter eode tssisg streaoaine ounrature aethod for the 
fl&it over ^ada to blade plai» of turbine casoades and subsequent 
applioatlofi of l^ls node for redesign of blade profiles * Tbe 
program is essentially f^ subs<»iic :^ow exeept for some local 

t 

reglmm of.supeisonic flow* trass osio flow solution is 

obtaiimd by a combination of a finite difference, stream function 
solution and a velocity gradient solution. For the solution of 
second degree equation of stream function the finite difference 
teehnique is used at the reduced weight flow. Subsequently the 
veloci^ ^^ient appros^h aalies use of the results of this analysis 
for full weight flow calculations. The program input ii^ludes 
cascade geometryt inlet and outlet flow angles and sta^iation 
flow conditions and provide the blade surface velocities and 
streasline oo-mrdinate in passage output. The computed blade 
surface velocities have been compared wi*^ the available experi** 
mental data. A systematic computational design system for 
ibirbimt cascade is i^^porated to re-sign and optimise the airfoil 
geometry to improve the blade surface velocity dislu^i^tios* 

The modified profiles thus obtai^d show improvement over the 
original ones and are libely to have lower profile losses due 
to reduced suction surface diffusion. 



CHi^TEE I 


INTROIXKITIOJi 

In 05 ^« 3 f to pi^dlct the pesrfo3rffience of & turbc^achine, m 
analysis of flow through it, is required. Tha internal flow 
through tiirhosiachines is usually three dimensional viscous 
and unsteady. The extreme coio^lexity of tne fi«H« mates it 
difficult to determine the velocity distribution and flow direc- 
tions throughout the flow field, simultaneously taking into 
account ail the factors, 'ihe turbine design oroblHu thtrefore is 
usually divided into three 12a jor steps, the first in the deter- 
mination of the requirement of tne flow, v^ork and speed. Too 
second step is the evolution of the velocity diagrams which should 
be consistent with the desired efficiency and nua&or of stages. 

The third step is the design of the blading that will provide 
the necessary flow angles and velocities required by the velocity 
diagrams* 

The essence of the classical approach is thh splitting up 
of the three dimensional probli^n into two separate two dimensional 
problems. The first problem is that of the flow In the aeriodional 
plsoie having r-z coordinates. The plane extends from hub to tip. 
The second plane is the cascade plane also known as blade to 
blade plane having the r-© coordinates. The Fig. l.l shows l^ie 
blade to blade and meridional planes as Si and S 2 surfaces 
respectively# 
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Through flow ec^ations have been deveiqped for these planes. 
Advent of a^dem ccffiqputers with their vast memory storage capa- 
bilities enables one to solve the extremely conplieated three 
dimensitmal non viscous equations. The computations using the 
cc^uters forms a framework for both rapid preliminary evaluations 
as well as for more detailed final designs* 


A successful dosigr for turbonochinexy cascades wuld 

provide a calculation tool b/ "Ahicn cascade geometry could be 
computed from certain design inputs and provide a subsequent 
procedure by which optimum geometries may be detexiained and chosen* 
The design procedure usually starts from known flow conditimas 
ahead of and behind ihe leading and trailing edges respectively* 

The thermodynamic states and flow velocities at the inlet and outlet 
of the cascade having been calculated in advance* The deirign 
prc^lem of a cascade is therefore to find out the blade profiles 
that lie cn a surface of revolution an.l also satisfy ii^osed 
design requirements* 


1*1 LITERATURE SURVEY 
1*2 IfCOIl^RESSXBLE FLOW 


For incojqpressibl© ficms, th© problmii is posed mathematically 
as the solution of th® Laplace equation 



in a potmtial field for specified boundary 


conditions* Direct problem approach for incoa^ressible potential 
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flow needed few specific Informations such as cascade geoi^tiryf 
inlet flow angles and discharge conditions* There were various 
methods available for the solution of the inconpressihle potential 
flows namely conformal transformation» singularity method etc* 

The conformal trmasformation method provides solution for isolated 
airfoil of arbitrary shape* It applied fCutta-JcmkimslEy transforma- 
tion for the solution near the pointed trailing edge$ and an integral 
equation was solved iteratively* This method finds notable appli- 
cations for cascades* The sljfigularity method had been develc^ed 
for thin airfoils with small carioer and pointed trailing edge* 
Theodorsen transformation has been used to determine the pressure 
coefficient around the airfoil* The resulting pressure distribution 
is corrected for interference effects due to other airfoils in 
the cascade* using suitidhle distribution of vortices and sources in 
in the singularity method* 

The pioneer of cascade theory was Chaplygin CD* A rapid 
appro ximatimi theory was given by Erwin and Yacobl (2) based 
upon known flew around an isolated aerofoil* ^ and Brown (3) 
gave a simple solution which is best suited to cl^ely-sbsoed 
tuidaine blading* Theodorsen (4) had produced a solution to ih9 
direct problem for isolated aerofoils of aj^ltrary shape* 
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1*3 COMPRESSISLE FLOW 

TUs inlet Mach Number for turbine blade rows usually has 
a high subsonic value and therefore a theoretical approach of 
subsonic coapressible flow through cascades is essential* Extensive 
work has been done for the case of compressible flows and various 
coHpressible flow theories have been developed* Matrix aethodt 
Finite element method^ Time marching method and streamline 
curvature lasthod have been used with success* 

1.4 MATRIX METHOD 

In the matrix through flow method the grid is fonaeil and 
the nonlinear stream function Is espressed in the form of quasi 
linear equations* The contribution made by Marsh (5) was in his 
selection of a very appropriate irregular finite difference grid* 
The finite difference equation is then formed and solved for stream 
function at all mesh points* The new stream function distrllmtion 
is liien used to form a new set of equations and the process is 
repeated until a convergence criterion is satisfied* As with ^e 
iterative procedures it is good practice to increase the inlet 
Mach Number in gradual st^s to the desired value* 

1.5 FINITE ELEMENT 4eTJK)D 

Finite elmaent ii^thod is more general aoid has advanta^ 
of flexibility as coiqp»ared to finite difference methods* It allews 
the use of m Irregular gild and higher order elements may be 
imposed locally to iiq>rove accuracy where boufidaries are irregular* 
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Finite element may be choosen ®u©h that grid points coincide with 
the boundarVf since the ccmptitaticn locations are also the comers 
of the elements ♦ The finite element method proceeds by ihe forma- 
tion of integrals of approximation over the domain under considera- 
tion* This is more readily achieved by identification of a varia- 
tional principle requiring maximization of a scalar functi<»i over 
the domain* 

One technique which seems to work for subsonic and low 
transonic flows is the local linearisation approach of Habashi (6) 

He linearise the equations with regard to the local Mach Huiai>ers 
throughout the region* A separate linearisation is therefore 
performed for each element t the variational principle for that 
element being based on the Laplace equation for the locally linearis^ 
cespressible flows* This linearisation is deployed within an 
iterative scheme and provides stsbility and rip id convergence 
for all subsonic Mach lumbers* Adler and ICris^rman (7) adopt 
a different approach by rewriting the governing equations in the 
form of Poisson equations having right hand side coi^osed of 
knoim functions* 

1*6 TIME MATCHING METIKID 

The reeimt develcpment in cascade flows in blade to blade 
plane includes time marching technique for the ihroughflow calcula- 
tions* For the case when there le locally siipersonic fl^i the 
region is no Imiger elliptic but is hyperbolic* The problem of 
Bix^» elliptic and hyperbolic domains with strong discontinuities 
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has been solved by use of techniques uslnq liie 

governing tiae dependent equations for coispiessible flow* This 
method demonstrated the ability to capture strong shocics* McDonald(8) 
assume isentropic flow in his finite volume tiro marching approach* 

Good agreement is obtained between isentropic time marching methc^i 
and e 3qoeria«>ntal data for cases having local Mach Number as high 
as 1*43* In an toproved scheme Denton (9) choose to solve -aie 
energy equations thus considering entropy vari^d&les in the flow 
field* Hovick (10) produced a predictor corrector time marching 
solution to the Havier Stokes equation and dei^^strated good agree- 
ment with eiqaerlmental results for a transonic flows* 

1*7 STREAMLINE CURVATURE METHOD 

The method is So nraed because the equatlcMss of motion are 
used in a form.t^lch relates the cross stream velocity gradient to 
the radius of curvature of the streamline* Since the radius of 
curvature of a streamline is a geometrical factor and stream 
spacing is determined by the c<mtinuity eqpati<m» it is relatively 
simple »thod to solve the equaticms of BK>tion for any giwn geometry f 
The advantage of these methods i^pear to be efficittscy of q^eration 
Mid the capability to consider velocitiew shove the critical* The 
methods have originated using streMnline and * quasi orthogonals* vhicif 
are simply lines passing from one channel wall to the other • in 
Ml arbitrary diiectlMi* It Is Mtvantageous to switch to full 
orthogcsial grid for transonic cascados* 
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In a typical atreaaliae ciii*^atura« i^og^amm tfae input 
data are differentiated twice and tde values obtained are 
used in Use velocity p*adient equation. A start is made witb 
approximate velocities along 'Uie aid streamline obtained frcmi tto 
previous iteration. A numerical integration is ^es performed in 
botb direotiom to find tbe blade 8urfa<» velocities, fhe velocity 
levels are ^len ad^Justed to satisfy continuity iterative ' teebnique. 
In i^rder to ensure rapid convergence an approximate damping factor 
must be obosen* Tbe ebasge in velooi’^ gradient ai^ars to be 
sensitive to tbe convergence criterion and tbis is reduced between 
iteratiotui to prevent divergent • 

Tbe streamline curvature technique is solved by two distinct 
teebniqii^s one using normals to the flow by Jansen (11} and 
Blnden and Carmichael (12) and the otl^r xuiing quasi-orthogonals 
by Theodore' Kats aula (13) and wilkiisson (14). The approach using 
normals to the flow has been xnsed tar two dimensional eases and 
has been less good ^lan qus»l*«or^ogonals In many respects. Quasi 
or^ogonals ate In practice near straight lines and proceed from 
CMse well to the other of ttm channel or one blade surfer to c^yaer, 
Gonfining Urn flow. All the mass now mosses the qussi-crl^ogonals 
and is used to satisfy the eontinulty equatloa* by ad;^stiiig tSm 
velocity along the quasi-orthogonals to obtained correct mass 
flow* fhi veloeily gradient relationship along the quasi ^thogonals 
is derived from momentum energy and entri^y equations. Since all 
the terns of right hand side are not known iterative techniques 
are imed to solve the slmaltaneous equations. This approach is 
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tmliCallj appa.ica1}I« to sulwofilG XLoms Inzt can also deal witb 
tr&EWOfiie flaws* 

Coi^sldaralJle work has heen done oa fiOQ>«eries airfoil 
dasipt over saries airfoil for highly loadad turbiot blade 
:i^ofiles by Konello and Mitohell (24), In t^ir work the non- 
series airfoils were designed from aeries airfoils and tl»y 
exhibit considerable improveit^nt over t\m aeries airfoils* 

2,Q. Ye (23) describes a ccmiiwtatiwjal design system for two 
dimension turbine cascade. He has made a sequence of calculations 
in idiich aii'foil profiles are designed frc® velocity diagram 
requirements and specified geometric parametei®, 

1.8 SCOPE (F PRESENT WMBflGATICm 

The present work incliided systematic computational design 
system for two diaei^ional turbii^ cascades using streamline 
curvature scheme for transonic flow* Considerable efforts have 
been made to implement the program given idlef* 21. As the program 
in Ref. 21 has been <teveloped after a series of modifications over 
the previoimi program given is Ref* 17t18,19 aikl 20, as such all 
the subroutii^s are not explained in a single source. Thei^ft^re to 
fully understand and is|4e»ent the program required considerable 
time and effort. The program has been run and tested for several 
cases, the computational results are compared wi-tti experimental 
results. There is good agreement between the results so obtained. 
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fhe work also laelodts a of oaloulations to modify 

to oxistifig Hade profiles os log the teobalqoe gl¥en ty Te to 

yield tetter Hade surface velocity distribution. Tbus in turn 
is reported tc reduce tbe pr-oflle losses over the blade surfaces 
and to eurb tbe tendency of early separation* For such investlga- 
tioiis three different types of stator and rotor blades were chosen 
and £&f each case velocity distributioii and pressure coefficient 
distribution over surfaces are calculated. Ihe results have been 
compared wiHi the results of the original profile. In all -Uiese 
cases the blades were highly loaded and the leading coefficient 
ranged from 0.8 to 0.94. The outlet MACH MUMSBR for rotor and 
stator blade profile ranged froo 0.^8 to 1.Q5. 
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CHAPTER 2 

FCBWHATIOH AHD APPROACH TO SGLCTI^ 


2.1 FORIRJLATIOii 

Tbe so3.utloQ fcs* tbe transooic flow is to lae obtained in two 
stages* In tbe fixst stage the oaXculatiens an done for reduced 
weight flow by finite difference solution for stream fuiKStion* The 
weight fl«nf is reduced sufficiently so ^lat the flow beccs^s subsonic 
-^oughcut the passage. 

In secoxxl stage the oaleulations are dcme with full weight 
flow using wlocity gsradient (streamline curvature)method* .The 
solution obtaii^d by reduced weight flow is used as input for 
full weight flow in velocity gradient solution* 

The basic reason behind tdie si^ittlng of the solution in two 
stages is the limitatifMis of tiie respeetive methods* The finite 
difference solution of the stream j^usstion equation is limited to 
strictly subsonic flmis* On the other hand a simple velocity 
gradient method is limited to a well guided ebannel* The purpose 
hezein is to combine these me-^iods so as to extend ^ reni^ 
of oases* 

The fcill owing simplifying assuaptioms are used in deriving 
am equaticms for both stages • ' 

1 • The fluid is perfect gas 

2* The flow is simady relative to the blade 

3* The flow is irrotational 
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4* fM HmMl is 0 (»'«viseou 8 and tisera is no beat transfer 
3» fba lilaSa to blade surfaoe is a surfa<» of roYolution 
6 * fbere sbonld not be my velocity oonponent normal to tbe t^ade 
to blaie surface* 

7« fbe stagnation temperature is uniform across tbe iMet 
8« tbe magnitude and direction of velocity across u|»treea 
and doenetream ^mddaries are to be uniform. 

9. the only forces are those due to masenttm and pressure gradient 

10. the relative velocity is subsonic f^ tbe first stage 
calodLations only* It may ise either subsonic or supersonic for 
the Mcond stage calculatic^# 


the soardiMte systne is shom in V Fig. 2.1 « the variables 
r and z are not totally independent on a stre^ siiirfa<^ and one 
of these variables can be elininated* therefore it is better 
to use neriidional coordinate *m* as indapendeot variable in place 
of r and z. tluis *b* and *9* become the two basic independent 
variable* Blade to blade surfaiMr of rev^ution is shown in Fig .2 *2 
and flow in nixed flew stream channel is show) in Fig. 2«3v 


For Umi nathematical f emulation of the problem Hut stream 
fheotion is used* tbe basic differential eepation which must be 
satisfied by t]t» atream function uad^ the gii»ii aisimptK^ is 
j ^ ^ 1 1 if ig (5lj»4.1 \ 

• 5 ? an ) an 


2 b 


siiK 


( 2 . 1 ) 
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fim desisriptiOB and derivation of the above equation is given 
in Ref. 15* 

fhe streamf unction u in to be normalised so that it has the 
value *0* on Itie upper surface of the lower blade and *1^ on ti^ 
liraer suz^ace of upper blade, fbe derivatives of the streamfunction 


satisfy 


( 2 . 2 ) 


au bf r 

V 


(2.3) 


The boundary conditions used for the rmgitm ABCIEFQR for 
finite flow region is shown in Fig. 2.4 are as folloi«. 


Boundary se^nt 


Boundary condition 

- u is 1 less tbe value of u on GH at 

the sase m oo-ordinate 

• tt *« 0 

-> u is 1 less than tbe value of u on £f 
at tbe sane m coordinate 


- (^ ) 

out 


out 
^ ^t 


• u is *1* greater the^ the value of u 
on CD at tdie same m eo*Qrdinate 

» u • 1 

- u is *1* greater ^aan Ibe vali;^ of u on AB 
at tbs sase *m* coordinate 

- (ia ) . . 

in ^ * ijB 



13 . 


Tb® equation (2.1) iJs no longer elliptical for locally 
superw^io jQ.©es but it becomes byperbolic. In such regions 
velocity graiieiit approach has been used for calculations. Simplifiaid 
velocity gradient equation is given below 


*a AW 4- B 


(2.4) 


where 


A * cos% + sinc^ tan |3 (1 cob^r ) 

^ dm 


(2.5) 


is used m blade surface 
2 ( ?Sbm 


sin 


( 


^u 

15 


m 

da 


AiE \ 2 
\Ldm I 


+ sioot. tanp (l+coB^ ) 

m d u 

J5 (2.6) 


is used at interior poiala md 


■m. ^ ^ dw 2 r gipio< 

B » r tanB ■■■■ —” ■ ■ + •.=— i-js=si^ 

° vttMp ^ T cosp 


(2.7) 


the description and derivational . of the above equations is 
given in Appendix *A* • fbe quantities used in above equati<»^ 
are kamm if solution of eqi^tion (2.1) is iinown. Tfaerefore, by 
using finite difference technique an apfroximate solutiom for 
reduced we i^t flow is obtained, such tiuit ^ streesaines f(sr the 
two flows wro nearly cl^e. For ^ faeoapressible flow there 
should not be any ohange in stresnlires if o) Is reduced in the sare 
prop<^tiOQ as w. This can be seen in equation (2.1) , since eoeffieient| 
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ar« eorntmt for inocitpressiliLe flow and equatioi} (2.1) does not 
ehac^ if is kapt constant, fbere will be change in strecatline 

simpes for coapressible flow case and factor which is propor- 
tional to weight flow for incoapressible ease is not sase strictly. 
HCwerer» for mmy cases tt^ error iatro^iced by using quantities 
in equation (2.5) ^ (2.6) and (2.7) fi’oia reduced weight flow 

solution is not significant. 

the eepaticn (2.4) can be solved nmerically along *0* 
60*#rdiiiateCi»e » at ocmstant m)* by using approxiaal^ value of A 
and B calculated by finite difference technique for completely 
subsonic flow. The continuity equation for flow through cascade is 

©2 

j f W coop b r d© » w (2.8) 

where is the value of © at the lower boundary and ©2 
is ^e value of © at tiie upper biamdary* 

Sinee tiMi unhnown variables in equation (2 .4) are not unique 
therefore it requires m iterative solution with the e(»idition 
it Bust satisfy ccntinuity equation (2*8). fo start the iteration 
sone initial value of W is taken in equation (2.8) and then 
eqimtion C2«>4) is solved* this procedure is related wi^ the 
new value of w until the correct sdutl^ is obtained. 

2*2 FIMHE BirSSiSiCB i^mMMkTim 

Finite difference apprcaination has been used for the solution 
of attem fcnction u* For s-hrean function solution a rectangular 
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grid bas astal^isfaed. Finite difierenee teohiiique has to be 
applied at each node points of rectangular iprid. Hon-linear stream 
fuBCtiGSi equation 




df 


du 




« £|£s?sifiiA 


can be written for each node points to Caloilate unkncnai stream 
fdncticm -yiroughout liie regi^. the ecpations are nonlinear since 
the coefficients depend upcm density. For n number of unknowns 
there are *n* ncmlinear equations. These nonlinear equatic»is may 
be sdlved iteratiwely by suc^ssive over-relaxation as described 
in Ref. 16* 


A typical mesh point with the msibering used to indicate the 
mesh poinlm is showa in Fig. 2.5. The equation (1) can be 
approximated by 


2u^ 2 u2 2U3 2u^ 2u^ 

iqtEpjT * ■ -E;^ * 



“2 ■ “1 


) * (-Siaio 


^4 f 


4 

? 



w 


b ? Sin 


(2.9) 


where 


h, . ( e), 

hj - r„ ( 9)2 
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ftoe eottf^deiits of io equation (2.9) ar® calculated as 


foUM. 


>^0 * ^ ®i ^*1 *’■ 

1x1 


.wiiem 


‘iz • • *34 “ * ®o ■ *12 * *34 


^ ^2 - Pl V \ ^4 * **3 SlnKo 

'>12 W¥'" ' "S'* *-TOW- --r^ 




< - '>12> > *2 


^12 ^ 

a " '~' ^1 

m 


j a^ 


^34 

*0 


and k- » 
o 


bo ^o 


sla «i 


( 2 . 10 ) 


Tl» equAtieeis eaa be used at all interior stesh points and 
for mAh poiAts ad^acMtat to Hade surfaces * Along tbe btouadary 
value of u is unlAom* the equation ulll vary* Alcmg upstream 
»«>u»»ary ) U too*, and . flnlf diff.r.iic. .pprxtaatlon gl«.. 


u^ *s u^ + h^ 


inlet 


inlet 


inlet 


( 2 . 11 ) 


and along the dovmstre^ boundary 


t. /»U \ „ »„ / *^^out 1 

- “3 * »3 ‘-35 - “3 - *>3 ^ 5 


(2.12) 
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The eal<ailati<»is of streaun function al^g -tiae houndary as shown 
in Fig. 2*6 can be derived by using periodic boundary conditions. 

Point *1* is suicide the boundary, Frosi the condition of periodicity 

where point IfS is above 1 at distance 6. Substituting this 

condition in equation (2*10) gives 

4 

Wq * SjjU^ g + *” ^0 

* i*2 

The value of streiei function u along MS. as shown in Fig. 2.7 
is *1* greater than the correspoiviiag valtMi along AS. Thus 

“2 ** ^ 2-8 ^ 

where point 2*s is distance *s* below MM in negative *8* 
direction. Therefore substitution in equation (2,10) gives 

“o • ®l‘l + '^“2-«.»3“3+®4"4 + ®2 ♦ '‘o <2.13) 

The D equations are reiresen'ted in the matrix form as 
Au ®* ^ (2.14) 

where ■ 

u » Ug, — Uq)^ 

where all eoBponeats are ualaiown values. 

A is the coefficient aatrlx of equatl^ (2,10) to (2,13) and 
whosa components are known constants. 
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Hit soiliU'tJ.cxi of oqiififtioti (2*10) by fteralrive procedure 
i« giveo fey 






a+1 


^+1 




(2.151 


for i s 1,Z, Q 


wfeert 


Is over*relaxation factor. 


The 


are tbe eleaeot 


of malariac A, and are eoapc^nts of Jc, 

fl» optiMia over'-reiaxatioa factor which aaxiaizes the asymptotic 
rate ©f coQvergtfice ©a© be exiaessed equivaXeotly as 


1 -f(0) 


fhe above equatlcm to determise optimum over’-relaxatloi} 
factor is derived ia ilef. 1i. 


2.3 CAucuuTioN Of mocm sm mmin 

Oace the viCLues of stream functiorm^ are ofetaioed oeiog fimite 
difference approxiaati^t the ^rivatives ^ and ^ calculated. 
The relative velocity W is calculated using relationship 

+ w| 

aM for each new value of W, W is calculated. The product If 
has its maximti® value W « relative velocity M has 

two values ©^ is subsonic and the other ©ne is super8onj.c. The 
subsonic valmes of relative velocity are ealeulaied corres ponding 
to the given viaues of w using Sewton’s method. 
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fbe tbtrffiodyiiaiBic relatioiahlps used io calculatioos are as 


1*- « 1 

IMMMMIIK a I 

in 




2 C fl 
p io 


irith the assuaptloo of iseo tropic flow 

f „ Vi> -1 

"j" '" iwrr* } 

y io ^in 


prodtiet P¥ is ealeuiated as 


i__ 
u -1 


the derivative of aliove relatioo wi-y^ respect to w is 

2»v 

.2 \ 5^ 




in 

S' in C’ - 


P "-ia 


4- 2'Ox- 

2^*in 


1/V-1 




The first value of relative velocity W is 


« 


tf^lgi 




S* islet 

The e<|uatioi} for liie aeart value of W hy liewtOB*8 aetiiod is 




I w • w 


a 


where a » 0,1 »2» 


( 2 . 1 « 


(2.11 


<2.li 


(2.1| 


(2,29 
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2*4 CALCUIATICK 

OF CRITICAL VELOCiry 

or 


J2. 2>)R 

■ ' 

T« 

(2. 2-1 

where t« * 

2uix» (to r)^ 



( 2 . 2 J 


2.5 CAUGUIATICN Oi* KAXIMOM VALUE W MASS FLOW W 

The st£PdJiisB vfdue of smss flow is attaioed wben 

p 2 i / V “I 

wj- + 2«A- (wr)'^ \ 

» ?’iii Wgj* (l 2 Cp / C2.23| 

2*6 mmiAfiQu m i^adihq miLm edge radii 

If m* and r* are tbe eo'-ordiiiatas of eeatre of tee are ^ben 
(m-ia»)^ + (© - 0*)^ » R^ (2.24 

wboro R is the raAiiis of leading or trailing edge. 

2 #7 CALOJLAIIOR OF SWFACl LESGIH 

fhe surface «l: leng12i of the hlade is calculated by the 
relaticsishiF 

n - . 

^ y ^i * ^®i * ^ f iz .2% 

2.8 VELOCIK CaaADIEHf APHlOACH 

Velocity gradients are the lines that extend fron one blade 
surfais^ to imother surface* The veloci% gradient equation is 

II - A¥#B 
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Aq IsitiAl tstlnate of W on the lower boundary is available 
frcm tbe reduced weight flow solutl^e. Tiae initial velocity 
is obtained by dividing the reduced tesigbt flow by reduced wei|^t 
flew factor* A nuaerical sduti^ to equation (2.4) is calculated 
by Ewiga«Ctitta apprtmb as ftOlowss 


.th 


If W|^ is known at i'” point calculated by 


^!^i « % ^ <^i«i <»i^i - ®i> 

bf . , + w*f. , 

Then «* 


The weight flow can be calculated by 

% 

f»l eo« p br <» 

The procedure la repeated until the calculated vNslght flow 
coitfieides with the actual weight flow. 
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2.9 AIRFOIL GEOMETRY 

In th® deslgpa of turbomachinery cascades, it is highly desirable 
that the curvature of the profile contour be continuous in order to 
avoid sever# velocity changes Ref. 17 and 18, In cascade design 
it should be possible to introduce local changes without substantiallv 
affecting the remaining part of the profile. The airfoil profile 
is fitted by segments of n^ degree polynomials while maintaining 
continuous curvature variation. 

The basic profile data for the design problem as shown in 
Fig. 2 . 8 . are suction and pressure surface defining 

points. A cartesl^ co-ordinate system is used with the Y-axis 
passing through the cascade leading edge and the x-axis through 
the centre of leading edge circle* 

Tht governing equations for determining the polynomials axe 
described' by I.Q.Ye in his paper Ref# 23* These equations are 
siodifitd according to the requirement of TSCWIC pTOgrat&me* The 
resulting equations ares 

^S1 - 

Bf 

** ^1 ** 2 ^ 

a*, 

^si * * "" 2 ^ 

Y 53 m s -Cs sin P cos ( r Bj) - H 

^3 ** 831 ^ (f+BJ) - R^, 

^S3 * ‘*‘^2 ** 2 ^ 

p# 

Yg^ « Sin (p* - H 
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0* 

%5 w [l-cos (P| - -^)3 

% « tan (P| - — - |) 

Pt 

Yp^ Bt •* sin (p£ +■ ' 2 "') 

* % [i+COS -f^)] 

^ P 1 “ ^"1 ** ^1 ’ ”* 2 '" "') 

^P4 * Sin (P| -f — } 

p* 

^4 * ^x"'^2 CP2 ^ ’ * *2 ' " ^ 3 

P* ■ 

^P4 t -—-2 *. J ) 

* ** <KC 4» ^ 

In the design polynoialals of degrees three and four are used 
to fit segment® of the profile to the set of ^^laetry condlticms 
deseriblng the oasoade* 

2*10 OETERiyill^ATION OF INPUT FCm THE 

The oatcade design can be Initiated by seleetlng a pitch 
to axial chord ratio S/p fro® loading considerati«Nns# The value 
of loading coefficient Is selected fro® design experience* 

From Horlock (25) the expreasion for is given below 

yj/j » 2 . Sin^ Pg (cot P^ +• cot Pg) 

Typical value for V'^ may rang# from 0*6 to 1.2 with the final 
choice being determined by the satisfactory performance of tho 
cascade* 
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A value of H can now be selected which deteiriBines the castbex* 
distribution and stagger of the airfoil. The cascade .gaging angle 
CPa ** ^ ('T/S)) is determined from the given exit gas angle ^2 

and nui^er The gaging length or cascade throat is 

calculated from 

A « S Sin 

The inlet metal angle coincides with the inlet gas angle at 
zero incidence. The exit iMtal angle and the uncovered turning P 
are finally decided on by finding an acceptable pressure distri- 
bution. 

The v.'edge angles and the auxiliary points (^ 2 * * 

"^54^* Tp 3^ ^ig* 2*8 on suction and pressure 

sides are determined in a trial and error procedure to give m 
acc^table airfoil shape* These variables strongly influence 
the cross channel area variation and h«nc# the pressure distributic«n. 



CHAPTER 3 


CCMilPiri ER PROGRAM 

3,1 INTRODUCTION 

Ihe second degree nonlinear equatlcms described In previous 
chapter are the governing equations for the blade to blade cascade 
flow and are solved by nusierical aslliods using ©odern coi^uters. 
Computer program for blade to blade flow have been develq[>«d by 
Katsanls in a series of publications Ref, 17,18,19 and 20, The 
latest version of Katsanis program is available in Ref, 21, 
published in 1%9* The present pr ogrma fonmilatlon is based on these 
publications of ^atsaiUs, The program is further suppleaiented by 
using plot-10 facility available In DBC-i090 system at X,I»T«ICanpur« 
It is used to plot velocities, blade shapes, surface pressure 
distribution and to cos^are and modify airfoils wife respect to 
existing airfoils* 

The computer program requires geometrical parameters as Input 
in mr# co-ordinates, appropriate gas constants mod operattog condi- 
tions such as inlet tesperature and density, inlet and i^tlet 
flow angles, weight flow and rotational speed. Fig. 3.1 m4 Fig, 3.2 
show the Oh# co-ordinate system for a typical bla#i to blade 
surface of revolution. The output obtained fiem ttte program intiudtos 

stream function values throughout the blade to blade passage, 
streamline locations, blade surface velocities and velocity magnitude 

and directions at all interior »sh points. 
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The whole program is to run in two stages. In the first stage 
the solution of nonlinear finite difference eguation is obtained 
through inner and outer iterations. Optiroum value of over^relaxation 
factor is calculated for Inner iteration to da^en the fluctuations 
and to lead the solution towards etmvergence# For each outer 
iteration the change in density is calculated and convergence 
of the solution is controlled by the variable Dentol {defined in 
input) • The solution obtained through first stage for reduced 
weight flow is used as input for s©C(^d sta^e to calculate ttit 
true velocities for actual weight flow with the help of velocity 
gradient technique* 

3*2 DESCRIPT 10^3 OF IHPUT DATA 

The variables used as input are geometric and nongeometxic* 

The gecesetricaX input variifeles me shown in Fig. 3*1 and Fig* 3*2 
The finite difference grid used In the program is showi in Fig*3«3* 

All the input varliybles are described below* 

AR * Gas constant J/tg/®lt 

BETAI -Inlet flow angis with respect to meriodional direction aloit 

BG 

BETAO- Outlet flow angle respect to rt^riodional direction 

along CF. 

BBTI1,BET12 - Angles with respect to merlodlonal direction at tangent 
points of leading and trailing edges radii with two blade 
surfaces » also tan P r (dO/dns) in degree. 
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BETOii BET02 - Angle® with z«®pect to i^rioslional direction at 
tangent point® of trailing edge radii in degree. 

BESP • Array of stream channel normal thickness corresponding to 
ME and ESSP arrays in meters* 

GHOEDF- Overall length in meriodit»ial direction, meters 
DBiTOL- Tolerance in d entity per iteration for re<feieed weight flow 
- Specific heat ratio 

MBI - Nuatoex of vertical mesh lines from AG to BG inclusive 

MBO - Number of vertical mesh lines from AH to CF inclusive 

MM • Total nuBfcer of vertical meshlines in omiri^ional direction 

from AH to BE. 

MR - Airay of merlodiimal *©* co*^rdinates of spline points for 

stream channel thicknessj meters 

MSPl, MSP2 Arrays of i!»riodi<wial ’m* co-o^wJinates of spline points 
on two blade surfaces, measured from the leading edge, meters* 

NBBI « Number of mesh spaces in ’8* <Slrectlon between AB and®. 

NBL » Nutaber of Blades 

NRSP - Number of spline points for stream chimmel radius (RWSP) and 
thickness (BESP) t co-ordinates* 

OMEGA- Rotational speed *w» rad/sec* 

ORF - Value of overrelaxation factor to be used in the solution 
of inner iteration 

3 ' ' 

RfK3IP- Inlet stagnation density, kg/m 



28 


REDF^- Factor by which weight flow wust be reciuoed in order to 
astore subsonic flow throughout the passage* 

RI1» RI2»Leadlng edge radii of the two blade surfaces, meters 
RG1,E02 “Treillng edge radii of the two blade surfaces, meters 
RMSP •• Array of ixo-ordinates of spline points for the stream 
channel radii correspcfiding to array, mtteri, 

STCBF - Angular ^-co-ordinate for centre of trailing edge circle cf 
blade with respect to the centre of leading edge circle of 
blade^radians 

SFLN01,SPU402 - MuaAjer of blade spline points given fur each surface 
as inputs include dui^es points at leading and trailing 
edge* 

TIP - Inlet staipation teraperature, 

THSPl,THSP2 - Arrays of ^8* co-ordinates of spline points correspimdiil 
to MSPl and MSP 2, r«iians* 

WTFL - Mass flcr# per blade for straa® channelf kg/sec. 

3,3 DE5CHIPT ION OF NORMAL OUTPUT 

The output Is controlled by the nusierlcal values of various 
variables used within the program* A value of *0* for any of 
these variables will cause the output associated with the corresptmd- 
Ing variables omitted. A value of U* will cause the corresponding 
output to be printed for the final Iteration* *2* for the first 
and final iterations; and »3Vfox all iterations. 
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AANDfC- Coefficient array constant vector and indexes of all 

adjacent points for each point in finite difference aesh* 

BLDAT— Ali geoiMStrical information which is constant 'Uiroughout 

the program (i.e* coordinates and first and second derivativei 
of all blade surface spiint points, blade coordinates » blade 
slopes, blade curvature, radii and stream channel thickness, 
blade surface angles and slopes, ITV and IV arrays). 

ERSOR - Maximum change in stream function at any point for each 
iteration* 

IHTVL • Velocity and flow angle at each interior mesh point for 
both reduced and actual weight flow* 

SLCEO - Stxeariiline ’S'* coordinates at each vertical mesh line 

STRFN - Value of stream function at each unknown mesh point in region 

SURVL « m co-^-ordinate, surface velocity, flow angle, dlstarjce along 
surface, and based on meriodionai velocity components, 

3,4 OVERALL PROGRAM LOGIC 

The compu ter program is composed of a main progrua and twenty 

seven subroutines, A list of these subroutines along with their 

functions Is as followst 

1, XRPUf - the subroutine reads all input data, calculates constants 

and iiltialires arrays, 

2, PHECAL • This subroutine calculates all quantities which remain 

constant for a single problem. 
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-• This ®\ibiroutin@ calculates the entries of the matrix A 
and the vector These coefficients must be recalculated 
for each outer iteration* 

4* SOR - It finds the linear solution to eciuation Au = ^ with fixed 
Coefficients by the method of successive over-relaxation* 
it also estimates the optiimifn over-relaxation factor* 

5* SLAX - It calculates the ctreamilne locations and / Wm 

6# TAMG - It calcuietes f amd then /V» and angle p throughout 

the region 

Im VELtXiy - It calculates the dimsity / and velocity % thr©u#iout 
the region and on the blade surfaces 


8* TVELCY - It solves the velocity gradient equatlran along each 
vertical line* 

It cem^^utes the coefficients acid constants along a 
vertical mesh line from blade to blade 
It calculates the values of h,r, and b. 

It computes coefficients and the const«fiis ki 
at a single point* 


9. COEFBB 

10* HRB 
11. AAK 


12* BDRY12 - Subroutine BimYi2 is called by COEFBB, It alters the 

values of h and r calculated by HRB for point 1 or 2 (Fig *11 


13. BDRY34 - Subroutine BDRY34 is called by COEFBB* It alters the 

values of h, r and b calculated by HR® for point 3 or 4 
{Fig.- 2.1)* 
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14* SLAV «• SLAV calculates and 3t7@ffifaliiie locatliHis* 

It also calculates at each aesh point and ^ on 
blad© surfaces* 

15* VEL Iha s^axliiaiM ralativa change In density alon^ a blade 
surface 1$ calculated in VIL* It also calculatts the 
value of density /, relative velocity W and ratio of 

16* SLCD « Subx^utine Bl^D calculates the ©-co-ordinate and d©/Cba 
of a blade surface for givcKii value of m, 

17* IPF - Function IPF nunbered the raesh point in the solution 
region* 

18* VELGRA - It calculates the values of constants A and B to solve 
the velocity gradient equation* 

19* MHORIZ - Subroutine MHORIZ calculates the n^o-ordinates of 

intersections of all hoxisontal raesh lines with a blade 
surface* . 

20. Da^SiTY- Subroutine PENSITY calculates the subsonic relative 
velocity W and corresponding density /, that result 
In a given value of the asass flow paraoieter fin* 

21* ROOT - Subroutine ROOT finds a root for fCx) « y by Newton’s 
isethod* 

22* SPLINE - It solves a tridiagonai aatrix equation to obtain the 

coefficients for the piecewise cubic polynostlal functlmi 
giving the spline fit curve* 
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23* SPL!^22 - It also solves the triadiagonal equation exc^t that 
for the end conditions, the slopes are specified, 

24# SPLIfIT - This subtoutine' is based on the cubic spline cumre 

The cubic spline curve is then used for interpolation* 

25* SEARCH -It calculates the mass flow parameter ^ W on blade 
surfaces# 

26* CCNTIN - CONTIH calculates a new estimate for the Initial 
value of W for equation C2*4), 

2'7# IMTGRL - IRTQRL calculator is the integral of a function passing 
through a given set of points. 

An overall flow diagram for tiie program is given in Fig, 3*4, 

The overall ccntrci of the calculation procecfeire is maintained by 

the main routine* 



CHAPTER 4 

RESULrS A'lD DlSCUSSlOfJ 

Th0 conputer progiraiB <31$euss®{| in the pre¥i<iHii chi^ter is 
applicable to axial, mixers ana radial flow turbines, since the 
effect of changes of cascade area in the flow direction and 
Changes of radius are included in the progxariu Howver, in the 
present appllcati.ons it has been applied for axial flow turbine 
cascade case only, including the rotor as well as stator cascadNs 
geometries. The test case have been chosen from Ref* 21 and 23, 
Initially considerable effort and time was devoted to obtain 
results similar to that of Ref* 22# This was necessary to validate 
the implementation of this program in the present investigations* 

The airfoil is redesigned by selecting airlbil cairt>er and 
thickness distributiisis # The ^plication of computer graphic 
design »thod# by 2,0* Ye (23) permits the designer to optimise 
the airfoil camber and thickness distribution, to minimize the suction 
surface rate of diffusion* The geometry redesign program is operated 
in an interactive mode by which the input geometric parameters are 
adjusted* To facilitate the initial specifications of a profile 
the few input parameters are Initially set to zero* The cubical 
polynomial subprogras SPilNT Is used to generate a pressure and 
suction surface blade profiles using the limited number of specified 
data* The profile is then successively modified by assigning values t<| 
the dummy points. Finally a profile is selected, after several 
iterations, after obtaining the ddslred blade surface velocity 
distributions. The velocity distributions Is normalised in terms 



34 * 


of critical velocity* The velocity distriiwtiofi so obtained has 
lower suction surface peak velocity and the rate of diffusion between 
the location of velocity peak and the trailing edge is lower ^an 
that in the case of original profile* The uniformity in velocity 
distribution is also obtained* 

TEST CASE I 

For the first test case ^e design parameters for stator 
guide vanes were taken fra® Theodore Katsanis {21)* The inlet 
and outlet flow angles assigned to stator guide vanes are *0* degree 
and *-67* degree respectively# The moderate flow has 0*:^ and 1*05 
M/yCH MUMBEii at inlet and outlet respectively* 

In the Fig* kJLA the analytical results obtained through 
program are cos^aredwith the ejgjerimantai results of J* Warren HWiitneyCsj 
There is a good agreement between the cc^uted critical velocity 
ratio with the experimental data* There is soim discrepancy near 
the trailing edge relion of the suction surface* The nunerical 
prediction methods are usually not so successful in predicting flow 
profiles in the vicinity of the blade leading and trailing edges# 

In Fig* 4*1B the stator guide vane profile from Ref* 21 
is shown by solid line* This profile is th^ modified by a systematic 
computational design system keeping all the thermodynamic parwters 
same* The modified profile is shovn by dashed line* There Is 
substantial ch«ige In th#rofll® shape of the suction surface near 



ih# txailing ®d 9 ®» Th© Fig* 4*2 displays th© surface velocity 
distributions* Th# isodified results have shov^ ir^roveiwBfit over 
those of original blade profile particularly near the trailing 
edge region# Further the laodified results are closer to the erperi** 
fflintal data* The pressure coeffici^t curves in Fig* 4*3 shows 
similar ia^roveaent* This iB?>roveraent over the blade surface is 
likely to fflinlmlze the profile losses ®od curb the tendwicy of 
separation on the suction surface* 

TEST CisySE 2 

In test cate 2 for stator guide van® th© geometrical input 
data are taken frm Ye (23)# The profile was assigned a loading 
coefficient of 0*8 with zero degree incidence and rnoderate subsonic 
exit flew (M^ ea 0*58, Wm.„*<D. 62)* The exit flow angle is -55®* 

In Pig* 4*4 the stator guide v«an® profile is shown by solid 
line# The jaccified profile Is shown by dashed line. It is displayed 
near th# leading edge region for both the pressure and cuciion 
surfaces# The coioparison of velocity ratio for both cases Is shewn 
in Fig# 4*5. Th© surface velocity variation in th® original case 
of Res, 23 are smoothed in the ©odified design of the present case* 
Similar type of ifflprovement seen in pressure coefficient distributions 
fox both surfaces as sho’/m In Fig# 4*6* Such an iisproved velocity 
distribution is likely to delay the seperation on the suctiork surface* 
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TEST CASE 3 

Unlike the first two i;ases the test case 3 is the turfsine 
rotor cascade redesign case. The rotor profile was assigned a 
loading coefficient of 0,94, The inlet and outlet flow lefigles are 
45.5® and -68® respectively. The moderate subsonic exit flow {^*.^=0,74), 

Fig* 4,7 is the con^iiuter generated graph of the turbine fcklade. 
There is substantially changes in camber and thickness distrlbiition 
in case of asodified profile particularly on the suction surface. 

There is a marginal reduction in the blade profile thickness in the 
modified case. The velocity dlstributi«n for both cases are showi 
in Fig, 4*8# The velocities obtained through modified profile for 
turbine blade show a good irijroveisent over prevlmis case. On the 
suction surface the peak velocity is lower In the modified design 
and algo more uniform velocity distribution Is obtained in the 
trailing edge region. Likewise the pressure coefficient distribution 
along the profiles are shown in Fig# 4.9, The modified results 
indicate the refinement of pressure coefficient distribution over 
previous OS'? and reduce the possibility of early separation because 
the minlaum pressure level is higher. 



4.1 CO!'iCLUSiaiS 


Th€^ first part of this work consisted of an impl omen tat Ion 
of the TSOhXC program for folarJe to blad.e calculations, giv^n in 
Hef, 21. Since this program in Hof, 21 has been devel^^<3 
successi've taodiflcations over previous programs given in Hof, ^7 3^3 j[_g 
and 2D* initially substas^ti&l tiiiie arid effort v'ss spent ij^ 
organising and Itrplementing it* The results presently obtained 
were checked wltn the results given in Ref. 21 for the sUndi^j.^ test 
case of s tatox blade row. After complete agreetwent was obtained 
in the two cases* otiwr test cates were taken up* 

In the second part of the work a systematic coaiputatlonal 
design systeia was used, rapid and effective d redesign of existing 
blade profilet turbine vascades. The approach follo^ved w®® 
similai' to that of Ve, The modified profiles thus obtaineci after 
several iterations show icsprovement over original one® In 
of bl'Ttde surface velocity distrli^^tion and are likely to encounter 
reduced profile losses. 

4.2 SCOPE OF FUTURE ViORK 

In the present vork efforts are made to modify the blade 
profile in order to reduce the profile losses. Du® to limitations 
of time however the effect of losses on surface velocity distil twtion 
hat not been included in the program. This can be done by * 

boundary layer analysis pxograpj. A program for reanalysis 'of leading 
and trailing edge velocity distribution can also be incotperated to 
obtain more accurate results near the loading wd trailing edge regiemf 
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DERIVATICW OF VELOCITY ORW3IEHT EQUATION 


Th. «locltwrx.<ll.„t ,<w.tl 0 „ i, „ e,p«,.lon of the force 

equ.tion. By e balance of force In the B-directlon, the foUowlng 

#qu«tion eaan bt obtained. 


H .1 »r^ 


(M) 


Tha tiiae darivatlve indlcataa the change in the quwitity for a 
moving particle as a function of time, m make use of the 
following relations (Fig* 2*1) 

- W Sin p 
r - W Cos p 
W » w sin a 

XU 

W « W cos a 
z la ' 

iSiC att W' 
dt r 

' dt " *6 

H ” *m 

H - » 

f -“f 

dt ■ ^S5 * "m S 
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W. c.n p.rfoi» th. InAcaUd dlff,rentl.tl«, of th. right side 
of .quotion (Al) by u,l„, th. pr,0Ml„, x.X,tl<«., «d then solve 
for aw/» (.*.lch eppeer. on both sld.. of th. «p..tlon) to obtUn 

^ dl * p) + r ton p ^ + 2w (A2) 


W. d..lx. now to .volnot. dp/dS to Uxm of first «.d SKond 
derlvotive. of 9 with «.p.ct to > .long streonllnos. Th. following 

relations hold for strsaalintt. 


tlUfi p « 



{A3) 

CA4) 


Equation CA4) is obtained by usinq equations (2# 2) and {2*3) 
Alio alonq streamlines we haw 

- Cos p (AS) 

** Sin tt (i^b) 


Now differentiate equation A3) and use equations (A5) m4 Ws) 


to obtain 



r Cos^p 




Along the surface of the blade d^/^^ cmi be ewily 


(A?) 


calculated since 9 is given txbliciliy as a function m* Heeieirer, 
in the passage d^/dm^ is glv«i Indirectly by the streaa function. 
Hence we will need an eipression for d%/d®^ in terms of the partial 
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strives of liio stroam function# PlxTit# firoro equation (A3) and 
(A4) wo havo 



(A8) 


By difforontiatin^ ocfuatlon (A8) wo obtain 



(du/d®) ^ CA9) 


Finally by using th# fact {fim oqs* CA3) and (A8)J that 



We can obtain 



By using equationa (A7) and (A30), equation (A2) can be put in the 


following foi» f 



« AWtB 


(All) 


whore 

A m cot^ «in a tan § (l+cos^ p) 

ctei^ 


(Ai2a) 


is used on blade aujrfaeof 


CM la,# 
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^ (A12b) 

If Dit«d at Intarior points md 

B ** r tan P (A135 

Equations (All) to (A13) ara in the foira used in th® prog^raa* 
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Fie 1.1 SCHEMATIC VIEW OF BLADE TO BLADE PLANE SI AND 


MERIDIONAL PUNE S2. 
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FIG 2.2 BLADE TO BLADE SURFACE OF REVOLUTION 
SHOWING M-0 CO-ORDINATES 











FIG 3.1 GEOMETRIC INPUT VARIABLES ON BLADE TO BLADE 
STREAM SURFACE. 
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FIG 3.2 GEOMETRIC IhPUT VARIABLES DESCRIBING STREAM CHANNEL 
IN MERIDIONAL PLANE. 
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FIG 3.4 FLOW CHART FOR SUBROUTINE CALLING SEQUENCE. 
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AXIAL POSITION X/BX 


FIG. 4,1 A COMPARISION OF BLADE SURFACE 

VELOCITIES WITH EXPERIMENTAL DATA 
FOR STATOR CASE 



THETA 



FIG. 4. IB COMPARISON OF BLADE PROFILE 










W/WCR 



AXIAL POSITION X/BX 
FIG A.S COMPARISON OF BLADE SURFACE 

VELOCITY FOR STATOR GUIDE VANE 







FIG 4.7 COMPARISION OF BLADE PROFILE 
FOR TURBINE ROTOR CASE 






PRESSURE COEFFICIENT 



